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Abstract 

We investigate inflationary models leading to density perturbations with a spectral index n > 1 ("blue 
spectra"). These perturbation spectra may be useful to simultaneously account for both the amount of 
ultra large-scale power required to fit cosmic microwave background anisotropies, such as those measured by 
COBE, and that required to give bulk motions and structures on the ~ 50 Mpc scale. 

Subject headings: 

Cosmology: theory - early universe -galaxies: formation - large-scale structure of the universe. 
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1. Introduction 



The recent detection by COBE (Smoot et al. 1992) of anisotropies in the Cosmic Microwave Background 
Radiation (CMBR) at large angular scales has provided an insight on the primordial density perturbations. 
This measurement makes it possible to normalize the power spectrum at the largest observable scales and 
constrains the primordial spectral index in the range n = l-lj^o.55- This is consistent with the well-known 
Harrison Zcldovich spectrum (n = 1), widely used in models of structure formation, but also with a sig- 
nificant range of values for the spectral index around it. It is however clear that, to further constrain the 
shape of the power spectrum, other measurements at smaller scales are required. One proposal (Davis et 
al. 1992) is to combine the COBE detection with measurements of the CMBR at the degree angular scale, 
as tested by the South Pole experiment (Gaier et al. 1992; Schuster et al. 1993). It has been claimed that 
these measurements are consistent with the COBE data for a scale-invariant power spectrum and the Cold 
Dark Matter (CDM) scenario; but the present accuracy and the inherent cosmic variance of the observations 
do not allow a precise determination of n (Dodelson & Jubas 1993; Bunn et al. 1993; see however Gorski, 
Stompor & Juszkiewicz 1993). 

On the other hand the standard n = 1 CDM model, normalized at ultra large scales to COBE (<~ 10 3 /i _1 
Mpc), predicts too much power on scales of the order of l/i _1 Mpc (where h is the Hubble constant in units 
of 100 km s" 1 Mpc" 1 ), which is in disagreement with observations such as the pairwise galaxy velocity 
dispersion (Davis & Peebles 1983). Tilted (i.e. n < 1) CDM models have been proposed to alleviate this 
conflict (Vittorio, Matarrese & Lucchin 1988; Cen et al. 1992; Lucchin, Matarrese & Mollerach 1992; Adams 
et al. 1992; Tormen et al. 1993). Other alternatives are based on either invoking a relic cosmological constant 
(Efstathiou, Sutherland & Maddox 1990; Turner 1991; Kofman, Gncdin & Bahcall 1993) or adding a hot 
dark matter component to the CDM scenario, leading to the so-called mixed dark matter scenarios (Davis, 
Summers & Schlegel 1992; Klypin et al. 1992). The latter possibility seems to be among the most promising 
alternatives to standard CDM, although at the expense of the introduction of one more free parameter in 
the theory: the relative amounts of the two dark matter components. 

Piran et al. (1993) have recently suggested an alternative view, which is to combine the CMBR mea- 
surement at large scales and the observation of large voids in the galaxy distribution. Since the voids can 
be assumed to result from underdensities in the true matter distribution, they can be used to normalize 
the power spectrum on the typical scale of the voids detected in the CfA survey (De Lapparent, Geller & 
Huchra 1986; Geller & Huchra 1989), i.e. of order 50/i _1 Mpc. This is independent of the galaxy to matter 
distribution relation, i.e. of bias. This analysis seems to favour a spectral index n ~ 1.25. An n > 1 ("blue") 
spectrum, normalized to COBE, might also provide the large power shown by galaxy peculiar velocities 
on similar scales, but it exacerbates the problem of the scale-invariant spectrum: the excess power on the 
Mpc scale. However, as pointed out by Piran et al. (1993) the extra power on small scales implies that 
non-linearity should be important on scales larger than usually thought, which must be properly taken into 
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account when studying this range of scales. Another possible way out of this small scale problem is to 
invoke a hot dark matter component, which would wash out the excess small scale power by free streaming. 
Moreover, the possible presence of a gravitational wave contribution to the COBE detection (Lucchin et al. 
1992; Davis et al. 1992; Lidsey & Coles 1992; Salopek 1992; Liddle & Lyth 1992, Souradeep & Sahni 1992) 
allows a further shift to the blue side of the spectrum, if the above w 50ft- -1 Mpc normalization is taken at its 
face value. This is because a stochastic background of gravitational waves only affects the CMBR isotropy, 
through the Sachs-Wolfe effect, on scales above w 1°, i.e. above 100ft -1 Mpc. 

On the other hand it is also possible to put upper bounds on n. A recent study (Carr and Lidsey 1993) 
has put a limit on the spectral index, n < 1.5, from primordial black hole overproduction constraints on 
nucleosynthesis and the COBE observations. 

Power spectra with n > 1 have not been very widely studied in the literature, one reason being that 
they are not the most common prediction of theories for the origin of the primordial fluctuations. Our aim 
in this paper is to explore inflationary models that originate this kind of power spectra. We first investigate 
in the next section models that can give rise to an exact power-law spectrum putting emphasis on the n > 1 
case, and then in section 3 some more realistic models where an approximate power-law with n > 1 arises. 
Finally we present our conclusions in section 4. 

2. Inflationary models with exact power-law spectra 

Quantum fluctuations of scalar fields during inflation are one of the most appealing explanations for 
the origin of the primordial density perturbations. An inflationary period in the early universe gives rise to 
both scalar and tensor perturbations in the metric. The scalar perturbations, associated to fluctuations in 
the energy density, are often described in terms of the gauge-invariant variable £. Its power spectrum to first 
order in the slow roll approximation is given by 



For the gravitational wave spectrum we have 

P Gw( k ) = ^ 



(1) 

aH=k 



(2) 

aH=k 



Chaotic inflationary models based on a polynomial potential give rise to a power spectrum P^ 2 {k) 
which decreases with k but only logarithmically. Another example is the power-law inflation, which is based 
on an exponential potential for the inflaton; it gives rise to a power spectrum decreasing with k as a power, 
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P^' (k) oc k a , with a < (Lucchin & Matarrese 1985). Power spectra decreasing with k are the most 
usual outcome of inflationary models involving one minimally coupled scalar field. This is determined by 
the dependence of H 2 /\<f>\ on a in equation (1). Or, equivalently, as by definition during inflation d increases 
with time (a > 0), we can say that d(H 2 /\<j>\)/dt < gives rise to a power spectrum decreasing with k and 
d(H 2 /\4>\) I dt > to one increasing with k. The fact that the time dependence oi H 2 is usually the dominant 



one and H 2 is always decreasing with time makes decrease with k. However, this is not a necessary 
conclusion: inflationary models in which increases with time more rapidly than H 2 decreases will 

originate a power spectrum Pq increasing with k. We want to analyse the possibility that a power spectrum 
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increasing with k as a power arises during the inflationary dynamics, that is P^' (k) oc k a with a > 0. This 
corresponds to a spectral index n = 2a + 1 > 1 for the density perturbations on constant time hypersurfaces 
before matter-radiation equality. 
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According to equation (1) a power spectrum P^' oc k a results when 

ff 4 



A(aH) 2a , (3) 



A2 



with A = H^J \4> 2 a 2a H 2a ) , where the subscript * indicates that the quantities are evaluated at the time a 
given wave- number fc* crossed the horizon fc* = a* if*. Using the equation of motion 2H = —K 2 <j) 2 , with 
k 2 = 8ir/Mp, we obtain 

Using the expansion parameter a as a "time" variable we can write equation (4) as 

where ' = d/da, which can be easily integrated. 

We first consider the case O^a^l. From equation (5) we get 

/ k 2 a- I x V(«-i) 



with 



In the a = 1 case we get 



and for a = 



13 ,-^StZ±^ +H ^. 











Va 2 


«§)) 



ff = ff exp — -j--, , (8) 



K 2 



H -- H ° + A 1 ^)- < 9 > 

We can now obtain 2 (a) by replacing H (a) from equations (6), (8) or (9) in equation (3). Using again 
a as variable we can integrate it to obtain 0(a) 

H 1 ~ a (a) 

d<t> = ±da ™ w . (10) 

For the three cases ^ a ^ 1,q = l,a = 0, the integration can be analytically performed, and the result 
inverted to obtain a (</>). With all these elements it is possible to determine the inflaton potential that 
originates these spectra. Replacing <p 2 (a) from equation (10) and H(a) from equations (6), (8) or (9) into 



we obtain V(a) and, as we also have a(<fr), we can obtain the inflaton potential V{4>). Since the computations 
depend on the value of a considered, they have to be analysed separately. 

For O^a^l the integration of equation (10) with if (a) given by equation (6) depends on the sign of 
(a — I) I a and B. We are going to analyse in detail the case (a — 1)/q < and then outline the main results 
for the other cases. The considered case corresponds to values of a in the range < a < 1, or to values of 
n in the range 1 < n < 3. From equation (7) we see that in this case B > 0. Integrating equation (10) and 
inverting the resulting expression one easily obtains 

„2a _ K 2 (l-a) 



2aAB 



■cos- 2 (V0, (12) 



with ip = (j)^J \n 2 \a{\ — a)\. Replacing in equation (11) equations (3), (6) and (12), the resulting potential 
is: 

V W) = "27i z^ 2{2 - a),{a - 1) {^){-a+('A-2a)^ 2 {^)) . (13) 

kt(1 — a) 

This potential has a minimum at ip = ir/2, two maxima, one to the right and the other to the left of the 
minimum, and it goes to — oo for ip equal to and n. The potential minimum is V m i n = 3B^=i / k 2 . (While 
writing this paper, we learnt that a similar solution has been obtained by Carr and Lidsey, 1993.) 

The solution obtained is inflationary, i.e. a > 0, only when <fi 2 < V{<p). Replacing equations (6) and 
(12) into equation (3) and comparing with equation (13) one can show that inflation occurs only for values 
of ip such that sin 2 (?/>) > a, that is in a region around the potential minimum compressed between the two 
maxima. Thus, inflation will occur while the inflaton field rolls down to its potential minimum, starting 
when its kinetic energy becomes small enough at values of 4> such that sin 2 (V') > a. From equation (12) we 
see that infinite e-foldings occur befor the inflaton arrives at the minimum: thus the field does not oscillate 
around the minimum and the model has no built-in reheating mechanism. As in the exponential potential 
inflation we have to invoke another process to end inflation. 

It is interesting to analyse the contribution of the gravitational wave modes in this model. According 
to Lucchin et al. (1992), the relative contribution of the tensor and scalar modes to the cosmic microwave 
background radiation anisotropies is given by 

(ST/T) S 2 H 



(5T/T) T 5 K 

In our case we obtain 

(ST/T) S 2 1- a 



(14) 

tan 2 (V>). (15) 



(5T/T)t 25 a 

We see that for values of <j> close to sin 2 (i/;) = a, where inflation starts, (8T/T)s/(ST/T)t ~ 2/25, thus the 
gravitational wave contribution can dominate; however, as <j> evolves towards the potential minimum, the 
contribution becomes always smaller. We can compute the number of e-foldings during which tensor modes 
dominate. From equation (12) one sees that the number of e-foldings it takes for <p to evolve from (f> x to <j> y 
is given by 

/ajA = 1 co^ 
\a x J 2a cos 2 {ipy) 
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Taking <p y as the value of <f> for which the scalar and tensor contributions are equal, we have from equation 
(15) cos 2 ip y — 2(1 — a) /(23a + 2), and <j> x as the starting value cos 2 tp x = 1 — a, we get 

.n(^=^.n(l + ^V (17) 



2a 

Thus, the gravitational wave contribution to the CMBR can only be dominant for the first few e-foldings of 
inflation. 

The other cases (a — l)/a > 0, a = and a = 1 can be solved following the same steps. We are going 
only to outline the main results. 

For (a — l)/a>0(a<0ora>l)we see from equation (9) that the sign of the constant B is no longer 
fixed and we have to consider the three possibilities: B = 0, B > and B < 0. The first case leads to an 



exponential potential V(cf>) cx exp(±K0^/2a/(a — 1)) and corresponds to the well-known power-law inflation 



case. Inflation can only happen for a < 0. The universe expands by \n(a y /a x ) = k\4> x — <p y \y/(a — 1) /2a 
c- foldings when the field evolves from <p x to <p y . The relative contribution of scalar and tensor modes to the 
CMBR is (5T/T) 2 S /(5T/T)^ = 2(a- l)/25a. We are now going to quote : a) the potential, b) the condition 
for inflation (a > 0), c) the number of e-foldings it takes the field to evolve between two values, d) the ratio 
of the scalar to tensor contribution to CMBR, and will discuss the main consequences for the other cases. 

Case B > 

1 

a) V(i/>) = ^Tq) cosh 2 ' 2 -"^"- 1 ^) (-a + (3 - 2a) cosh 2 ^)) (18a) 

b) tanh 2 (V>) < (a- l)/a-> a > (186) 

2/ 



1 / sinh (ip x 



{5T/T)% = 2 a-1^ 
{5T/Tf T 25 a 



The shape of V(tp) depends on the value of a considered. For a < (n < 1) it has a maximum at 
4> = and decreases monotonically on both sides, going to for </> — > ±oo. The expansion of the universe is 
inflationary for all the </> values; for large (f> the potential looks as an exponential, as in power-law inflation. 

In the a > 1 case there are three subcases. For a < 1.5 the potential has a minimum at (j> = 0, two 
maxima, one at each side, and it goes to — oo for <p — * ±oo. For 1.5 < a < 3 the minimum is at (f> = 0, there 
are no maxima, and V((p) — ► +oo for <fr — * ±oo. Finally for a > 3, = is a maximum, there is a minimum 
at each side and V{4>) — > +oo for — > ±oo. The inflation condition a > indicates that inflation happens 
only while 4> evolves towards <p = for \<p\ smaller than a given value and that it takes infinite e-foldings to 
reach = 0. Finally we see that the gravitational waves contribution to the CMBR may dominate, but only 
for values of <p close to the one in which inflation starts; this can last for few e-foldings, as equation (17) 
holds also in this case. 
Case B < 

a) V(i/>) = ^(^"i) sinh2(2 ~ a)/(a ~ 1) W (-<* + ( 2a - 3) sinh 2 (^)) (19a) 



b) tanh 2 (» > a/ (a - 1) -» d > (196) 
lnf^=^ln^^ (19c) 



a x / 2a \cosh (i/) y ) 

Condition b) indicates that there is no inflationary solution for a > 1. In the case a < the potential 
goes to —oo for \<j>\ — > 0, has two positive maxima and goes to for |<^| — > oo. Inflation only happens for 
large values of while it evolves to oo. In the large-0 limit this case reduces to the exponential potential 
case. 

The two remaining cases, a = 1 and a = 0, are solved in the same way, using equations (8) and (9) 
instead of equation (6) . 

a = 1 case 



/2 

b) |0| < ► a > (206) 

c) ln teH(f:) (20c) 

pr/r)| _ 4 i 

dj (5T/T)2 " 25 W [Md) 

The potential has a minimum at <p — 0, a maximum at each side and it goes to — oo for \<j>\ — > oo. Inflation 
happens while rolls to the minimum in a region between the two maxima. The gravitational wave contri- 
bution is dominant only for the largest values of \<j)\ satisfying the condition for inflation, and for In (25/2)/2 
c-foldings. 

a = case 

The case of an exact scale-invariant spectrum is given by: 

b) k\4>\ > 4V2 -> a > (216) 
ln(j)=g(^-^) (21c) 

d) (^VT)|_^ (21d) 
j (<5T/T)2 " 400 • 1 j 

It corresponds to the so-called Intermediate Inflation (Barrow and Liddlc 1993). The potential goes to — oo 
for \(j)\ — > 0, has two positive maxima and goes to for \(j)\ — > oo. Inflation occurs while <j) g° es — > ±°° f° r 
large |0| values. The gravitational wave contribution to the CMBR can dominate only for the smaller values 
of \4>\ compatible with inflation and for 23/4 e-foldings. 



3. Working models 

The exact solutions with a > presented in the previous section illustrate the fact that it is possible for a 
minimally coupled inflaton field to produce density perturbations with spectral index n > 1, provided that its 
potential has a particular shape. The main characteristic is that the potential minimum has a non- vanishing 
energy density, which corresponds to an effective cosmological constant. Some simple potentials with this 
characteristic originate (non-monotonically) increasing power spectra. For example V{(f>) = Vo + f<j> n gives 
rise to a spectrum increasing with k for wavelengths that left the Hubble radius when f<f> n < Vo, and 
decreasing with k in the opposite case. 

An obvious problem with these inflationary scenarios is how to produce the reheating. There are two 
main ways in which reheating is achieved in inflationary models: through rapid oscillations of the inflaton 
field around the potential minimum or through a first-order phase transition with bubble production. The 
first one is associated for example with a polynomial potential inflation model, while the second one is 
associated with some extended inflation models. In the cases discussed above, where an inflaton field rolls to 
a potential minimum with large cosmological constant, none of them works. Due to the large friction term 
in the Klein-Gordon equation, the scalar field does not oscillate around the minimum but approaches it 
with an ever-decreasing velocity, taking infinite e-foldings to reach it; thus the oscillation mechanism cannot 
operate. The first-order phase transition can work, in principle, provided that the cosmological constant 
corresponds to the false vacuum energy density of a second scalar field, which at a given time makes the 
phase transition to its true vacuum by bubble nucleation. However the problem with this picture is that, 
as the inflationary expansion of the universe is close to exponential, the graceful exit problem of the old 
inflationary scenario appears (a slower expansion as in extended inflation would avoid it). A possible way 
out of it (Linde 1991a, Adams and Freese 1991) may be to introduce a coupling between the fields, such 
that the tunnelling probability is determined by the dynamics of the other field. This reheating problem can 
nevertheless be avoided in a new class of models, the "hybrid inflation" ones, recently proposed by Linde 
(1991b, 1993). In these models inflation ends in the following way: the slow roll of the inflaton field at a 
given time triggers a second-order phase transition of a second field, whose false vacuum energy density is 
responsible for the cosmological constant term. The field rolls down to its potential minimum and oscillates 
around it. A particular potential in which this scenario can be realized is (Linde I99Ib,1993) 1 

VOM = ^(M 2 - Aa 2 ) 2 + ^ + ^V. (22) 

For <j> > <fic = M/g, a is in its false vacuum at a = 0. The potential for <f> in this regime is given by 
V{4>) — M 4 /4A + m 2 </> 2 /2. When the cosmological constant term dominates, this gives rise to an inflationary 
phase with a power spectrum P^(k) increasing with k. When <j> < (j> c , a rolls down to its minimum and the 

1 The fact that n > 1 in this model was also noted by Liddle and Lyth 1993. After completing this work 
a revised version of Linde 1993 appeared where this issue is also discussed. The results are fully consistent. 
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universe reheats. The constraints put on the model parameters in order that it works are discussed by Linde 
(1993). In order to have inflation for 4> > 4>c one needs M 2 ^> 2m 2 X/g 2 , and in order that inflation finishes 
soon after <fi reaches <f> c , one needs M 6 <C 3A 2 Mptn 2 /2ir 2 . We can now look at the resulting spectral index 
for the perturbations a = d In P^ 2 /d In k. 

The Klein-Gordon equation for the inflaton can be written as 

<A" + 3<A'+^^> = 0, (23) 

where from now on a prime denotes derivative with respect to In a. For 4Am 2 / k 2 M 4 <C 1, the solution can 
be written as 



f) ■ ,24) 

The power spectrum P^ 2 cx H/<j>', and it is a good approximation to take dink = dlnaH ~ dlna. Thus 

/ H\' H' ft 1 4Am 2 
^("Vj =H-ft^^M- (25) 

The spectral index is positive, as anticipated. To be consistent with the approximations done, this result 
holds only for small values of a. This presents no problems for values of n close to 1, such as those required 
to explain the large voids observed (n ~ 1.25, Piran et al. 1993), which correspond to a ~ 0.12, consistent 
with the allowed range of parameters. This value fixes the relation among m 2 and M 4 through equation 
(25). 

The amplitude of the perturbations that left the horizon around 60 e-foldings before the end of inflation 
(the ones relevant for structure formation) is given by 



1/2,, , V2^gM 5 ( AXm 2 



(26) 



(The last factor, which was neglected in Linde 1993, has to be included when we consider values of 
Ato 2 /k 2 M 4 <~ a <^ 1, but not <C 60). This amplitude has to be <~ 3 x 10 -5 in order to agree with the 
COBE results. This, combined with equation (25), finally gives the values of m and M: 

M^^10- 3 M P , 



m ~ J-TT 2 X/g 2 lQ- b Mp. 



For values of the coupling constants A and g of order 1 these values are consistent with the above 

>i/ 

C 



1/2 

constraints. Thus, this provides a working model in which a power spectrum P> increasing with k can 



arise and a successful reheating mechanism can be implemented to end inflation. 
4. Discussion and conclusions 
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The solutions obtained in section 2 illustrate the fact that it is possible for a single minimally coupled 
inflaton field to produce density perturbations with a spectrum k a for any value of a = const. The decreasing 
power spectrum (a < 0) is the known output of power-law inflation. Usually, spectra with non-constant 
negative a are the results of the inflationary potentials considered in the literature. The increasing spectra 
(a > 0) are less frequent. They arise in models in which \<fi\ decreases more rapidly than H 2 . This condition 
is realized in the cases discussed while the inflaton rolls down to a potential minimum with a non-vanishing 
potential energy (i.e. a cosmological constant-like piece). Mathematically it is difficult to write down a 
general condition for the potential to ensure that the power spectrum increases with k, but it is possible to 
obtain the condition in the slow-rolling approximation (<j> <C 3Hcj), and <j> 2 /2 <C V), that is VV^/V^ > 
3/2. An inflaton field with a potential of this type will generically have problems to reheat the universe. As 
we have discussed in the last section, a possible way to circumvent this is given by the "hybrid inflation" 
models, where the coupling of the inflaton with a second field at a given time triggers a second-order phase 
transition of this field, which reheats the universe. 

Another interesting issue to discuss is the contribution of the gravitational waves to the CMBR. In all 
the solutions with a = const> the tensor mode can dominate over the scalar one only for a few e-foldings as 
the relative contribution always decreases with time. The reason is easy to understand, looking at equations 
(1) and (17). To have a spectrum increasing with k, H/\<p\ must increase with time. As H decreases with 
time, H/\4>\ necessarily increases with time and thus, according to equation (17), the ratio of the scalar to 
tensor mode contributions increases with time. We have seen that the maximum number of e-foldings, since 
the beginning of inflation, during which the gravitational waves can dominate is 

1 / 23 \ 

N GW = -m(l + -<>). (27) 

Thus, if gravitational waves contribute significantly to originate the CMBR anisotropics at large angular 
scales, the corresponding wavelength should have left the Hubble radius not very far after the beginning 
of inflation. This corresponds to a choice of the end of the inflationary expansion in a specific range of 
values, which makes the hypothesis unattractive. The more general prediction of these models is a negligible 
gravitational wave contribution to AT/T. 

The results presented in this work show that power spectra with n > 1 can be easily originated in 
inflationary models. The astrophysical interest in these models comes from the necessity of simultaneously 
accounting for both the right amount of ultra large-scale power required to fit CMBR anisotropics and 
that required to give rise to bulk motions and structures on <~ 50 h^ 1 Mpc scale. A wider analysis of the 
astrophysical implications of this kind of spectra and a comparison with the available observations are needed 
to test the model further. 
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